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While for synchronous deterministic cellular automata there is an accepted definition of reversibility, 
the situation is less clear for asynchronous cellular automata. We first discuss a few possibilities and 
then investigate what we call phase space invertible asynchronous cellular automata in more detail. 
We will show that for each Turing machine there is such a cellular automaton simulating it, and that 
it is decidable whether an asynchronous cellular automaton has this property or not, even in higher 
dimensions. 

1 Introduction 

For synchronous deterministic cellular automata the topic of reversibility has gained a good deal of 
attention. Reversibility is known to be decidable in one dimension [ 1 ] and undecidable in two and more 
dimensions [4J. Reversible synchronous cellular automata are also computationally universal (see, for 
example, for the one-dimensional case). 

In the present paper we take a first look at the analogous questions for asynchronous cellular au- 
tomata. Hence the rest of this paper is organized as follows: In section [2] we introduce some notation 
used in this paper and two variants of asynchronicity. In section [3] we discuss several possibilities for 
the definition of "reversibility" of asynchronous cellular automata. In section [4] we show that each Tur- 
ing machine can be simulated by a phase space invertible purely asynchronous cellular automaton, and 
in section [5] that the property of being phase space invertible is decidable for several variants of asyn- 
chronous cellular automata. 

This paper is based on the diploma thesis of the first author ifTTTl . 

2 Basics 

2.1 General Notation 

We write N for the set of natural numbers without 0, Z for the set of integers, G Z d for the J-tuple with 
each component equal to zero, B A for the set of all total functions from A to B, and 2 M for the powerset 
of a set M. The cardinality of a set M is denoted as \M\. The restriction of a function / : A — > M to a 
subset BCAof its domain is written f\g. 

In this paper we are interested in <f -dimensional cellular automata 

(d £ N) and the set of cells is usually denoted as R, that is, R = Z d . If the set of states of one 
cell is denoted as Q, the set of all (global) configurations is Q R . A neighborhood is a finite set N = 
{n\ , . . . , nt} C Z d of (i-tuples of integers. A 
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local configuration is a mapping £: N — » Q; thus Q 1 ^ is the set of all local configurations. The local 
configuration Cf+jy observed by cell i G /? in the global configuration c G Q R is defined as 

c,-+iv :N Q, 

n i—T- c(i + n). 

The behavior of each single cell of a deterministic cellular automaton is described by the local transition 
function 8 : Q N —> Q. 



2.2 Asynchronous Updating Schemes 

A local structure (R, N, Q, 8) of a deterministic cellular automaton together with a prescription how cells 
are updated induces a global transition relation & C Q R x Q R describing the possible global steps which 
satisfies 

c^c' =>\li£R:c(i) G {8{c i+N ),c(i)}. 
We will use the same symbol for the related global transition function 

ch> {c | c&c'}. 

With this notation c S% c' is equivalent to c' G &(c) and both indicate that it is possible to reach global 
configuration c' G Q R in one step from global configuration c G Q R . 

In a global step each cell has two possibilities: to be active and make a state transition according to 
the rule or to be passive and maintain its state. Restrictions made by different updating schemes lead to 
different possible behaviors, that is, different relations/functions of cellular automata. 

Now, we will have a look at two different types of asynchronous updating. 



Purely Asynchronous Updating. This version of asynchronous updating has been considered for 
many years now Hl|3l. In order to distinguish it from the other form mentioned below we call it purely 
asynchronous updating. 

In each global step there are no restrictions on whether a cell may be active or passive. Thus in each 
step there is a subset A C R of active cells which make a transition, while the cells in the complement 
R \ A are passive and simply maintain their state. 

Note that A is is allowed to be empty. The additional requirement A / might look irrelevant, but a 
closer look at the constructions and theorems reveals that it would render the main results wrong (see for 
example the remark after Lemma [6]>. 

Given a deterministic cellular automaton and a set A C R we define the function 

Aa:Q R 
VieR: A A (c)(/ 

Synchronous updating is described by Ar. The union of all Aa, interpreted as relations, is the general 
step relation for purely asynchronous updating, for which we will write 




A= |J A A 

ACS 
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Note that for each global configuration c G Q holds 

A(c) = {A A (c) | AC*}. 

Fully Asynchronous Updating. In the fully asynchronous updating scheme it is required that in each 
global step exactly one cell is active. Using the notation from above one may say that one only looks at 
the relations Am where the set of active cells is a singleton. For the union of these relations we will write 

Ai = U A W ■ 

ieR 

Even for relatively simple deterministic cellular automata, for example, the elementary deterministic 
cellular automata or two-dimensional minority, the analysis of their behavior under fully asynchronous 
updating is surprisingly "non-simple" EH21H1. 

To distinguish the global step relations of two purely or fully asynchronous cellular automata C and 
G we will use A A , A, and Ai for the respective relations of C and Fa, T, and F\ for the respective relations 
of G. 

3 Which Definition of Invertibility? 

For the global transition function A of a synchronous deterministic cellular automaton C the following 
conditions are equivalent: 

Rl Each global configuration c G Q R has exactly one predecessor under A. 
R2 Each global configuration c G Q R has at most one predecessor under A. 

R3 The inverse of the transition graph of C, that is, the direction of each transition is reversed, is the 

transition graph of a(nother) synchronous deterministic cellular automaton C . 
For asynchronous cellular automata these are really different conditions. In that case, obviously, condi- 
tion Rl still implies R2 as well as R3. But the reverse implications do not hold. In order to clarify this, 
we first show that there is essentially only one asynchronous cellular automaton satisfying Rl or R2: The 
identity. 

Lemma 1. If a purely or fully asynchronous cellular automaton C = (R,N,Q,8) has a local transition 
function 8 which is non-trivial in the sense that 8(£) ^ £(0) holds for at least one local configuration I, 
then there are two different global configurations c and c and two singleton sets of active cells {a} and 
{a} such that Ar^i (c) = A{^ j (c). (Without loss of generality, we assume that is in the neighborhood.) 

In other words, there is a global configuration which has two different predecessors. 

Lemma[T]shows that the requirement of Rl or R2 leaves the trivial identity as the only asynchronous 
cellular automaton. Later, in section [4] and section [6j we will see, that there are non-trivial asynchronous 
cellular automata satisfying R3. Therefore, for asynchronous cellular automata, R3 does not imply Rl 
or R2. 

Proof (of Lemma^. Let £ G be a local configuration with q' = 8(£) / ^(0) = q. 

Let a G R be a cell large enough such that the neighborhoods (— a) +N and a +N are disjoint. We 
consider a global configuration c G Q R in which cell —a and cell a both observe £ in their neighborhoods, 
that is, C(_ a ) +A ? = I and c a +N = £, in particular c(—a) = q and c(a) = q. Define two global configurations 
c and c which are identical to c with the only exceptions c(— a) = q' and c(a) = q' respectively. Since 
q' 7^ q the global configurations c and c are different. 

But the two global configurations c' = A{ a j(c) and c' = A_ Q (c) are the same (see Figureflj). 
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Figure 1: Displayed are the relevant parts of the global configurations c, c, c, c' and c' in the one- 
dimensional case, that is, R = Z, with neighborhood N = { — 1,0,1}. Grey colored cells are in state q and 
black colored cells in state q' . The state of white colored cells is not further specified. The solid directed 
lines indicate the two interesting global transitions. 



• Each cell i ^ {—a, a} is passive; therefore c'(i) = c(i) = c(i) = c(i) = c'(i). 

• <*'(-£?) = c(-a) = q' = 8(1) = 8(c^ a)+N ) = c'(-a). 

• c'(a) = 8(c a+N ) = 8(1) = 4 = c(a) = c'(a). ■ 

In a short presentation given at Automata 2011 Sarkar and Das ifTUl also consider some kind of "re- 
versibility of 1 -dimensional asynchronous cellular automata". We note that at least their setting is 
completely different from ours: They only look at finite configurations (with both, periodic and null, 
boundary conditions). And as far as we understand their Definitions 1 and 2 iflOl p. 32], they call an 
asynchronous cellular automaton, restricted to configurations of a fixed length, reversible if each (finite) 
configuration has at least one predecessor. We do not pursue this line of thought. 

Instead we use R3 as the guiding light. In order not to overload the word "reversible" with too many 
meanings we will speak of phase space invertibility, which we will often abbreviate as invertibility. 

Definition 2. A purely asynchronous cellular automaton C = (R,N,Q,8) is called phase space invert- 
ible, if there is a purely asynchronous cellular automaton G = (R,M,Q,y) such that for each pair of 
global configurations c,c' € Q R holds 

c'eA(c) cer(c'), 

which is equivalent to 

3A C R : c' = A A {c) 3A' C R : c = T A ,{c). 
We give a name to the set of cells on which two global configurations differ in 
Definition 3. The difference D c c i of two global configurations c,c' £ Q R is defined as 

D cc , = {ieR\c i ^c' i }. 

Note that D c c i = D c i c . 

Since for purely asynchronous cellular automata there are no restrictions on the set of active cells, cells 
which do not change their state when active, can be removed from the activity set without changing the 
outcome of a transition. 

This suggests the characterization of invertibility given in 
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Lemma 4. Two purely asynchronous cellular automata C = (R,N, Q,8) and G = (R,M, Q, 7) are inverse 
to each other if, and only if, for each pair of global configurations c,c' G Q R with D c c i 7^ holds 

c = a Dcc , (c) c = r Dcc , (c). 

Proof. 

"only if": Let C and G be inverse to each other. Moreover let c,d G Q R be an arbitrary pair of global 
configurations with D c c j ^ 0. 

If c' = Ap , (c), then there exists a set of active cells A' such that c = r^(c') by premise. In this 
case D c c / C A', since c; = c- for each cell i G" A', and c = To ,(c'), since c,- = c- for each cell 
i G A' \D C . c '. Therefore 

c' = a Dcc , ( c ) => c = r Dcc , (c'). 

The other direction follows by symmetry, 
"if": For each pair of global configurations c,c' G Q R with D cc i / let 

c = a Dcc , (c) c = r Dcc , (c). 

Moreover let c,c' G 2 R be an arbitrary pair of global configurations. 

If there exists a set of active cells A such that d = Aa(c), then c' = Ad , (c), since D c c i C A and 
c- = c ; - for each cell i G A \D C C /. In this case we get c = To , (c r ) by premise if D c c < / or by 
definition of and equality c = c' if D c c / = 0. Therefore 

3AC/? :c ' = A i (c)^ C = r Dcc , (c'). 

The other direction follows by symmetry. ■ 

Analogously to definition [2] we define invertibility for fully asynchronous cellular automata in 

Definition 5. A fully asynchronous cellular automaton C = (R, N, Q, 8) is called phase space invertible, 
if there is a fully asynchronous cellular automaton G = (R,M,Q,y) such that for each pair of global 
configurations c,c' G Q R holds 

c'gAi(c) cGri(c'), 

which is equivalent to 

3a G R : c' = Ai a y(c) <^=^ 3a' G R : c = r^jfc'). 
We characterize inversion for fully asynchronous cellular automata in 

Lemma 6. Two fully asynchronous cellular automata C = (R,N, Q,8) and G = (R,M, Q, 7) are inverse 
to each other if, and only if, for each pair of global configurations c,c' G Q R with D C ( j = {a} holds 

C ' =&{a}{c) <^ C = F {a} (c') 

and for each global configuration c G Q R holds 

3a G R : c = A{ a j(c) <^> 3a' £ R : c = T{ a '} (c). 

A comparison with the formulation in lemma [4] shows a complication. This is due to the fact that for 
fully asynchronous cellular automata one always needs at least one active cell. 
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Proof (of lemma^ty. Let c,d G Q R be an arbitrary pair of global configurations. There are three cases to 
consider: 

Case 1: \D cc i\ > 2, that is, c and d differ in more than one cell. Then d ^ Ai(c) and c £ T\(c'), and 
hence 

c'eAi(c) <=> cenfc'). 

Case 2: D cc / = {a}, that is, c and c' differ in exactly one cell a. To reach d in one step from c by C or c 
in one step from d by G cell a must be active, thus 

d G Ai (c) <J=^ c' = A{ a } (c) and 
ceTi(c') <=^- c = F {a} (c'). 

Case 3: D cc / = 0, that is, c = d. If c' can be reached from c by C, that is, c' G Ai (c), c may not be 
reachable from d by G or be reachable with a different active cell, and vice versa. Thus, a better 



characterization of this case is not as simple as in the other cases and is postponed until lemma 15 



Note, that for purely asynchronous cellular automata this case is trivial: Simply choose the empty 
activity set. ■ 

Lemma [4] is only correct because we allow the activity set to be empty. If we would not allow this, we 
would additionally have to require for each global configuration c G Q R that 

3a G R : c = A{ a j (c) <^=^ 3a G R : c = Ts a n(c). 

As we will see, this makes it much more difficult to prove completeness and decidability. 



4 Turing Completeness 

We will show that invertible purely asynchronous cellular automata are computationally universal. 

It is known that reversible deterministic synchronous cellular automata are computationally universal. 
For one-dimensional cellular automata this can be shown by reversibly simulating reversible Turing 
machines which are computationally universal; for <f -dimensional cellular automata the result holds too 
(more details can be found in [5 ]). 

Hence any construction which transforms a reversible synchronous cellular automaton into an invert- 
ible asynchronous cellular automaton is sufficient to show the computational universality of the latter 
class. It turns out that Nakamura's method (H of transforming the local transition function of any syn- 
chronous cellular automaton into one for an asynchronous cellular automaton while "basically preserving 
its global behavior" (irrespective of reversibility) is all that is needed for purely asynchronous updating. 

In general the following transformations which maintain local synchronicity and guarantee invert- 
ibility are performed: 

1 . Each cell, additionally to its current state, remembers its previous state and manages a three- valued 
time stamp. 

2. Each active cell only changes its state, if thereby no information that may be needed by neighboring 
cells is lost. This is the case, if neighboring cells have the same time stamp or are one step ahead. 

3. Moreover, each active cell maintains its state, if the local configuration it observes is "illegal". 
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More specifically, let C = (R,N, Q,8) and G = (R,N, Q, y) be two synchronous cellular automata with 
corresponding global transition functions A and T. The interesting case will be that they are inverse to 
each other, that is, A -1 =F. Without loss of generality, we assume that the neighborhood includes and 
is symmetric, that h, N = {—n \ n <E N}. 

We will now construct an asynchronous cellular automaton C = (R,N,Q,8) from C. In order to 
save parentheses we will occasionally write £ n instead of £(n) (V/ £ Q 1 ^ ,n £ N) and c,- instead of c(i) 
(Vc £ Q R ,i £ R). The cellular automaton C is defined as follows: 

1. The set of states is Q = Q x Q x {0, 1,2}. For q = (q\,q2,t) £ Q we denote the first component as 
curr{q), the second component as old(q) and the third as time(q). 

2. Given a local configuration i we say, that cell is ahead (of its neighbors) if, and only if, there is 
an 7i £ ,/V such that time(lo) = time(£ n ) + 1 (mod 3). 

3. If in a local configuration I of C cell is not ahead define the corresponding current local C- 
configuration curr(£) of C (not C !) as 



Analogously, apply the following construction to G, resulting in G = (R,N, Q, f). 

1. Q = Q x Q x {0, 1,2}, curr(q), old(q) and time{q) are defined as above. 

2. Given a local configuration £ we say, that cell is behind (of its neighbors) if, and only if, there is 
an « £ A^ such that time(£~o) = time(£ n ) — 1 (mod 3). 

3. If in a local configuration £ of C cell is not behind define the corresponding old local C- 
configuration old(£) of G (not G !) as 



4. A local configuration £ is backward movable if cell is not behind and curr(£o) = S(old(£)). 

5. The local transition function S is then defined by 




curr(£ n ) if time(£„) = time(£o), 
old(£ n ) if time(£ n ) = time(£o) + 1. 





old(£ n ) if time (£ n ) = time (£q ) , 
curr(J n ) if time(£ n ) = time{£o) — 1. 




( old(£o), y{old{£)), time(£~o) — 1 ) if € is backward movable, 
In otherwise. 



We now have 



Theorem l.IfC and G are synchronous cellular automata which are inverse to each other, then C and 
G are purely asynchronous inverses of each other. 
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Proof. Let c,c' G Q R be two arbitrary global configurations and 

D = D c . c , = {ieR\ c t / c'i) 
their difference. According to lemma [4] it is sufficient to prove 

c'=A D (c) <S=^ c = T D (c'). 

We will prove in detail that 

c = A D (c) =>- c = T D (c'). 

Because of the symmetry of the constructions is it not surprising that a proof of the inverse implication 
can be given analogously. 

Now let c' = Ap(c). Consider an arbitrary cell i G D. Since c, / c • = 5(c I+ ^) the local configuration 
c i+ N is forward movable, because otherwise cell i would maintain its state by definition of 8. Therefore 

\/n £ N: time(a +n ) G {time (ci), time {ci) + 1}, 

old(ci) = y(curr(c i+N )), 
curr{c[) = 8{curr{ci + N)), 

old^c'j) = currifi), and 
time[c^) = timeifi) + 1- 

Now consider an arbitrary neighbor n £ N. Since we have assumed that ,/V is symmetric, cell i is a 
neighbor of cell i + n. 

There are two possible cases: 

Case 1: time{ci +n ) = timeici). 

Case 1.1: i+n G D: Then c ;+ „ ^ c' i+n and therefore 

time(c' i+n ) = time{ci +n ) + 1 = time{ci) + 1 = time(c'i) and 
curr(c i+N )(n) = curr(c i+n ) = old(c' i+n ) = old(c' i+N )(n). 

Case 1.2: i + n D: Then c, +n = c' i+n and therefore 

time{c' i+n ) = time{a +n ) = time(ci) = time(c'i) — 1 and 
curr(c i+N )(n) = curr(c i+n ) =curr(c' i+n ) = old(c' i+N )(n). 

Case 2: time(ci +n ) = time(ci) + 1. Since / is a neighbor of i + n and time(ci) = time{ci +n ) — 1 we have 
c' i+n = c i+n . Therefore 

time(c' i+n ) = time(ci +n ) = time(ci) + 1 = timeic'^) and 
curr(c i+N )(n) =old(c i+n ) =old(c' i+n ) = old(c' i+N )(n). 

Taken together we always have 

time{c' i+n ) G {time (c-), time (c-) — 1} and 
curr{c i+N )(n) = old(c' i+N )(n), 
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and therefore cwrr(c-) = 8(curr(ci + N)) = 8(old(c' i+N )). As a consequence c' i+N is backward movable 
and hence 

curr{ci) = old(c'j) = curr(Tr){c')i), 

old(ci) = y(curr(c i+N j) = y{old(c' i+N )) = old(T D (c')i), and 
time{ci) = time(c'j) — 1 = time(Fr)(c ),•), 

which is a long-winded way of saying c\ = To{c')i. ■ 

The above proof is incorrect if we restrict purely asynchronous cellular automata to non-empty sets of 
active cells. The problematic case happens when the minimal difference D is empty. In this case we 
cannot use it as the set of active cells. But nevertheless c' = c may be reachable from c by C in one step 
with a non-empty set of active cells and we need to prove that the same holds for G. We do not know 
whether this is always the case. 

One part of the problem is that even if c is a global configuration in which the registers old, current 
and time do not conform we need to show the property. This could be solved by adding more restrictions 
that do not corrupt the simulation as in the second part of the definitions of backward and forward 
movable. 



5 Decidability 

We will show that phase space invertibility is decidable for arbitrary-dimensional purely asynchronous 
cellular automata and one-dimensional fully asynchronous cellular automata by presenting two algo- 
rithms which always terminate and for any given automaton find an inverse if one exists. 

For any given asynchronous cellular automaton C these algorithms only look for inverses among the 



finitely many automata with the same neighborhood as C. This is justified in section 5.1 where we prove 
that if C is invertible there is an inverse with the same neighborhood. 

To decide whether two asynchronous cellular automata C and G are inverse to each other the algo- 
rithms only verify that C and G are inverse to each other on a subspace of the phase space that is restricted 



to global transitions in which only cells of a fixed finite subset of cells are active. In sections 5.2 and 5.3 
we prove that this is sufficient for arbitrary-dimensional purely and one-dimensional fully asynchronous 
cellular automata respectively. 



5.1 Inverse Neighborhood and Translation Invariance 

Consider any set of all global configurations that agree on cell and all its neighbors. For each of these 
global configurations make the global transition step where only cell is active. All resulting global 
configurations again agree on cell and all its neighbors. 

An inverse can undo all these transitions, whereby at most cell is active. Regardless of the states 
of cells besides cell and its neighborhood, the inverse does the same and can thus not benefit from a 
larger neighborhood. 

Definition 8. We say that a cellular automaton C = (R,N, Q,8) has minimal neighborhood if it has no 
dummy neighbors, meaning for each neighbor n £ N there exist local configurations £, £' € with 

%{»} = and ** + £ n SUCn that 5 W + 8 ( i ' )■ 

With this term we can rigorously state and prove 
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Lemma 9. Two inverse purely or fully asynchronous cellular automata C = (R,N,Q,8) and (7= (R,M,Q,j) 
with minimal neighborhoods have the same neighborhood. 

Proof (by contradiction). Assume N ^ M. Without loss of generality, let N\M ^ 0. Then there exists 
a neighbor n € N\M. Because the neighborhoods are minimal there also exist local configurations 
1,1 e (f with i\ N \ [h} = £\ N \ {A] and t fl / l h such that 8(t) ^ 8(1). 
Choose global configurations c,c £ Q R with 

q = 1% as well as Cf l = 

c\n\{h} = ^N\{n} = $\n\{A} = ^|jV\{n}» and 
c\r\n = C\R\N- 

For the global configurations d = A{ } (<?) and d = Ar } (c) we obtain 

d = 8(!)^8(i)=d , 

dfi 4 / £fi ^ni 
d\R\{0,n} = <^|ff\{0,n}> and 

d\R\{o} = c\r\{o} as well as <^Ir\{0} = c\r\{o}- 
Because C and G are inverse to each other we furthermore have 



c = Tn .(d) as well as c = Tn .(d). 

c.d v ' 

Note that D, t,D~ ? C {0}. See Figure|2jfor a graphical representation of the situation. 



A {0 } A {0} 



Figure 2: The labels on the dashed lines indicate on which cells the connected global configurations 
differ, the solid lines show possible global transitions. 



Case 1: n^O. Then co = Co. With do / do it follows that do ^ Co or do / Co- Without loss of generality, 
let do / Co. Then c = T^} (d) and thus co = 7(<^o+m)- 

Consider global configuration d £ Q R with ci|#\{o} = <^U\{o} an d do = <^o- From n^Mwe have 
^o+m = ^o+m a nd therefore 

T{o}(d)o = j(do+M) = y(do+M) = co = co and 
r {0(Wls\{o} =d\ R \{o} =d\ R \{ Q } =c\r\{o}, 
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Figure 3: The labels on the dashed lines indicate on which cells the connected global configurations 
differ, the solid lines show possible global transitions. 



which is a long-winded way of saying c = F{ Q } (d). 

From do = do ^ do = Co we conclude d = Ar } (c) = d and do = do = do, which contradicts do^do; 
see Figure [3] 

Case 2: n = 0. Then f. M and hence do+M = do+M as well as co+m = cq+m- 

Case 2.1: cq ^ do or <?o 7^ Jo- Without loss of generality, let Co / <io- Then c = r{ j(<i) and 
therefore co = /(^o+m)- With ^o+a/ = do+M it follows that <?o = /(^o+m) and therefore c 7^ 
c = r{ }(J). Hence c = Tq>{d) = d and c = r{ j(c). 

Finally, from cq 7^ co we have c = A{ j(c) and conclude d = d, which contradicts do / do', 
Figure [4] 



see 



{0} 



{0} 



{0} 



d\ 7^0 




Figure 4: The labels on the dashed lines indicate on which cells the connected global configurations dif- 
fer, the solid lines show possible global transitions, and the rounded rectangle signifies that the enclosed 
global configurations are found to be identical during the proof. 

,«. v V 

Case 2.2: Co = do and co = do. Then c = d and c = d. Consider global configuration b = T^j (c) = 
Tjo} (c). From cq 7^ cq it follows that bo 7^ co or bo^co- 

Without loss of generality, let bo 7^ co- Then c = A{ }(&) and with c/c also c = A@(^) = b. 
Therefore c = A{ }(c) = d = c, which contradicts do 7^ cq; see Figure [5] 

Because every possible case led to a contradiction, the assumption must be false, meaning N = M. ■ 

Since dummy neighbors can be added and removed without affecting the phase space we get 
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Figure 5: The labels on the dashed lines indicate on which cells the connected global configurations 
differ, the solid lines show possible global transitions. 

Corollary 10. For each invertible asynchronous cellular automaton exists an inverse with the same 
neighborhood. ■ 

Let us now briefly consider translation invariance. 

Definition 11. For each translation vector j G R the map Zj : Q R — )• Q R with 

V/ £R: Xj(c)i = c i+j 

is called (j '- translation. 

Similar to synchronous, asynchronous cellular automata are also translation invariant - we only need to 
translate the set of active cells as well. 

Lemma 12. For each translation vector j G R and set of active cells AC/? holds 

where A — j = {a — j \ a G A}. ■ 
5.2 Purely Asynchronous Cellular Automata 

If every transition of a purely asynchronous cellular automaton, in which only cell is active as well as 
arbitrary neighbors of cell 0, can be inverted by another automaton and vice versa, then these automata 
are inverse to each other. 

This follows from locality properties and translation invariance of purely asynchronous cellular au- 
tomata and is rigorously stated and proved in 

Lemma 13. Two purely asynchronous cellular automata C = (R,N,Q,8) and G = (R,N,Q,y) are in- 
verse to each other if, and only if, for each pair of global configurations c, c' G Q R with G D c c i C {0} UN 
holds 

c = A Dcc , (c) c = T Dcc , (c). 
Proof. The forward direction follows directly from lemma [4] For the backward direction consider any 
global configurations c,c' G Q R . 

First, let c' = Ad cc , (c)- Choose any cell i G D c c i that changes its state during the transition from c to 
c'. Because C and G are translation invariant it follows that 

r Dccl (c') i = T i (r Dcc ,(c')) 
= r Dc ^ i (z i (c')) 
= r Dcc ,_ ; (T I -(A Dc/ (c))) 
= r Dcc ,_ ! (A Dcc ,_ I (Ti(c))) . 
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With G D c l j — i we further get 

r Dcc ,(c , ) ! - = r {0} (Az )ff; - I -(T«(c)))o. 

Because the transition state of cell only depends on its neighborhood we farther deduce 

d 

r D cc ,( c ')i = r {0}{A(D cc ,-i)n({0}UN){?i(c)))0- 

V v ' 

d> 

We can make cells, which do not change state if active, inactive without disturbing the transition from d 
to d'\ 

d' = A Ddd ,(d). 

The difference D^j' of d and d' is contained in {0} UN. From G {D c c > — i) n ({0} UN) we furthermore 
get 

do = T,(c) = Ci ^ c- = 8(c i+N ) = 8(Tj(c) 0+N ) = d' , 
meaning G P>d,d'- Thus, by the premise of the backward direction, 

and therefore 

r o cc , ( c ')i = r {o} (d')o = r D( , d , (d')o = d = T(( c )o = 
Because cell i was arbitrarily chosen, c = Tp , (c') and we have proved the implication 

c = A Dcc , ( c ) c = (c). 
By symmetry we conclude that the implication 

c' = A D( . ( ,(c) ^c = r D(V (c') 

holds as well. Because the global configurations c and c' were arbitrarily chosen C and G are inverse to 
each other according to lemma |4] ■ 

Thus only cell 0, the neighbors of cell 0, and the neighbors of each neighbor of cell have to be consid- 
ered when testing whether two purely asynchronous cellular automata are inverse to each other, because 
the transition state of cell or an arbitrary neighbor of cell only depends in its observed local configu- 
ration. 

Theorem 14. Phase space invertibility is decidablefor purely asynchronous cellular automata. 

Proof. Let C = (R,N,Q,8) be an arbitrary purely asynchronous cellular automaton. Consider each of 
the 121'^'" purely asynchronous cellular automata G = (R,N,Q,y) in turn. Test whether for each pair 
c,c' G q{°} unu ( n + n ) of finite configurations with G D C( j C {0} U N the equivalence 

c = A Dcc , (c) c = T Dcc , (c) 

holds. If this is the case for one automaton G, then C and G are inverse to each other, and C is invertible. 
Otherwise C is not invertible. 

Note that N+N = {m + n \ m,n G A^}. ■ 
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Remark. It is in fact not necessary to consider all purely asynchronous cellular automata with the same 
neighborhood. One can show that if C is invertible, the local transition function of its inverse can be 
constructed easily from the one of C. Thus only one candidate needs to be considered. 
For further details see ifTTl section 3.5]. 

The above proofs are incorrect if we restrict purely asynchronous cellular automata to non-empty sets of 
active cells. But at least in the one-dimensional case invertibility is still decidable (see iPTTl section 6.2.3] 
for further details). This is due to the fact that a purely asynchronous cellular automaton that is restricted 
to non-empty sets of active cells is invertible if, and only if, it is invertible without the restriction and 
invertible as fully asynchronous cellular automaton (see IfTTl corollary 3.15]). And, as we shall see, the 
latter is decidable in the one-dimensional case. 



5.3 Fully Asynchronous Cellular Automata 

In each transition of fully asynchronous cellular automata exactly one cell is active, whereas in each 
transition of purely asynchronous cellular automata arbitrary cells are active. One may naively think 
that deciding invertibility is thus easier for fully than for purely asynchronous cellular automata. But the 
added difficulty is that in every transition one cell must be active. 

To decide whether two fully asynchronous cellular automata C and G are inverse to each other or not, 
we have to show that if C transits a global configuration c into d when cell a is active, that G can transit 
d into c when a cell a' is active or that no such cell exists, and vice versa. 

If c a 7^ c' a then a and a' must be the same cell or C and G are not inverse. But if c a = d a then a and 
a 1 are in general different cells. So how or where can we either find such a cell d or decide that no such 
cell exists? 

We answer this question in the one-dimensional case in 

Lemma 15. Let C = (Z,iV, Q, 8) and G = (Z,N, Q, y) be two fully asynchronous one-dimensional cel- 
lular automata. Define the maximal distance of a neighbor 



m ■ 



max„ e A?|?i| ifAf/0, 
iftf = 0, 



and a finite set of active candidate cells 



^ = |{-iei 2m+1 ,...,-i,o,i,...,iGi 2m+1 } ifw^-oo, 

[ {0} if m = — oo. 

Then C and G are inverse to each other if, and only if, for each pair of global configurations c,d G <2 Z 
with D c . d = {0} 

c ' = A {o}(c) c = r {0} (c) (i) 

and for each pair of global configurations c,c' € <2 Z with c = A^ j (c) and c' = Tjo} (c') 

3a G £/ : c = T {a} (c) and 3d G £/ : c = A {a , } (d). (2) 

In the proof of this lemma we need segments of global configurations: For each global configuration 
c G <2 Z and cells i, j G Z with i < j the -segment (of c) c[i,j] is the (J — i + l)-tuple with c[i,j] = 

5 • • • ,Cj). 
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Proof. If there is only one state, meaning \Q\ = 1, there is nothing to show. We thus assume \Q\ > 2. 

"=>•" Let C and G be inverse to each other. According to lemma [6} Equality ([T} holds. To show 
equality ([2]) we choose an arbitrary global configuration c E <2 Z with c = Arm (c) or c = r^} (c). 
At first, let N = ®. Then ^ = {0}. Let c = A{ j(c). Assume c / r{ j(c). Consider the global 
configuration d E <2 Z with c • = Co for each cell £ € Z. Then c' = A/m (c') and from 

c' a = c + r {0 }(c)() = T{ a }{c')a 

we obtain c' ^ Ts a \ (d) for each cell a E Z - in contradiction to the premise that C and G are 
inverse to each other. Thus our assumption is false, meaning c = T^{c). One can analogously 
show that c = T^{c) implies c = Ar i(c). 

Now let N y^®. The idea is to construct a global configuration d E <2 Z for which 

c 'eAi(d) <=> c'ET^c') 

does no? hold if there does not exist a cell a E such that c = r\ a j(c), if c = A{ j(c), or c = 
A{ fl }(c), if c = r^o}( c )' contradicting that C and G are inverse to each other. 

Consider the |g| 2m+1 + 1 segments 

c[a — m,a + m] 

for cells a E {0, 1, ... , |<2| 2m+1 }. Because at most \Q\ 2m+l of these segments are different, there 
exist two cells a, a' E {0, 1, . . . , |g| 2m+1 } with a < a' such that 

c[a— m,a+m] = c[a! — m, a' + m\. 

Analogously there exist two cells /3',/3 E {-\Q\ 2m+ \. . . , 1,0} with j8' < j3 such that 

c[P' -m,p' + m] =c[f}-m,P+m). 

P' P a a a' 

n I I I I I I I I I I I I I I I I I YA I I I I I ■ I I I I V/////A I I I I I I I I I I I I I I I I I 



Define for each index k E N 



and 



Ok = a' + (k — \){a' — a) — m, 
a' k = a' + k(a' — a)+m, 



f3 k = P'-(k-l)(P-P')+m, 
# = j3'-*03-/J')-i». 



II I II I I I I II I'll II III V//X/XA I I I I I I I I I I I I I I I I I 



p[,..., pi c «+" 
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Because < a < a' and j3' < j3 < we have 



and 



—m < a - m = oci < a 2 < a 3 < • • • , 
m < a' + m < a' + (a' - a) + m = a[ < a' 2 < otg < • • • , 



< A? < /3 2 < J3i = /3' + m < m, 

< $ < $ < J8( = P' - Q3 - J8') - m < J8' - in < -m. 



Consider the global configuration c £ Q with 



c'[j3' -m,a' + m] =c[fi' -m,a' + m], 

c'[a k , a' k ] = c[a - m, a' + m] (Vfc e N), and 
</\PM=c\P'-m,P+m]Qtk€N). 



, a?, a, 

a . 2 

Y//MVA I I MW/d I I MW/I I I I 



■K ft ft' ft c "+ n «!,...,«[ 



This global configuration exists, because the rear and front parts in which the definitions of two 
consecutive indices overlap are precisely the identical segments c[a — m, a + m] = c[a' — m, a' + 
m] and c[j3' — m, ft' +m] = c[/3 — m,/3 +m]. In fact, for each index k G N it holds that 

{j3'-m, ...,a'+m}n{ai,...,ai} = {a' -m,...,a' +m}, 

{p'-m,...,a'+m}n{p[,...,fr} = {p'-m,...,p'+m}, 

{a k ,...,a' k }n{a k+ i,...,a' k+l } = {a k+ i,...,a' k } 

= {a' k -2m,...,a' k } 
= {a k+u ...,a k+i +2m}, 

{/5;,...,ft}n{/5; +1 ,...,ft +1 } = {^,...,ft +1 } 

= {$,...,& + 2m} 

= {j3jt+i-2fn,...,j3fe+i}, 

and 

c'[a' — m, a' + m] =c[a'-m,a' + m] 

= c[a-m,a + m] =c'[ai,ai +2m], 
c'[/3'-m,/3' + m] = c[/3'-m,j3' + m] 

= c[j3-m,j3+m] =c'[j3i-2m,j3i], 
c'[a^-2m,a^] =c[a'-m,a' + m] 

= c[a— m,a+m] = c'[a k+ i,OC k+ \ +2m], 
c'[fi[,l5l+2m]=c[P'-m,P' + m] 

= c[P-m,P+m] =c'[p k+ i-2m,p k+ i]. 
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We will now show that for the global configuration d any cell whose neighborhood is not contained 
in srf behaves the same as at least one cell whose neighborhood is contained in srf if active. Note 
fhatAf C {— m, ...,m}. 

Assertion. For each cell a' G Z exists a cell a G {/?', . . . , a'} such that 

c'[a — m,a' +m] = c[a — m,a + m]. 

Proof. Let a' G Z be an arbitrary cell. If a' G {/?', . . . , a'}, choose a = a'. From now on let 
a' £ {/?', . . . , a'}. We consider only the case that a' > a! + 1. The case a! < j3' — 1 can be shown 
analogously. Thus let d > a' + 1. 

Because a' > oti + m and each of the intersections {at, . . . , a(} n {dk+i , . . . , } (\/k G N) and 
{a' — m, . . . , a + m} contains 2m + 1 elements, there exists an index k G N such that 

{a'-m, ...,a' + m} C {a k , . . . , a' k }. 

Choose a = a' - k(a' - a) G {a, . . . , a'}. Then 

c [a — m,a' + m\ = c[a — m,a + m\. u 

With this setup we can conclude the proof: 

1. Let c = A{ j(c). From 

c'[—m,m] = c /;;./;; 

we have d = A{ j (c ; ). Assume that for each cell a G stf it holds that c ^ (c). Choose an 
arbitrary cell a' G Z. Then there is a cell a G {/3', . . . , a'} such that 

c'[d — m,a +m] = c[a — m,a + m]. 

Because cell a is contained in it holds that 

c'a> =Ca+ F{a}{c)a = T {a , } {c') a > 

and thus c' ^ Y^ a i^{c'). Because a! G Z was arbitrarily chosen, this contradicts the premise 
that C and G are inverse to each other. Thus our assumption is false, meaning there exists a 
cell a G £/ with c = F^ a j (c). 

2. Let c = rr }(c). Like above one shows that there exists a cell a's^/ with c = A{ a /j(c). 

Therefore equality (|2]) holds as well. 

For each pair of global configurations c,c' G Q R with D c c < = {a} holds 

c' = A w (c) T (c / ) = T «( A w( c ))= A {()}(T a (c)) 

e l=f T fl (c) = r {0} (T a (c')) = T fl (r w ( c ')) 
c = r w (c'), 

and for each global configuration c G 2 s holds 

3a G : c = A{ a } (c) 
=>- T«(c) = T fl (A {fl} (c)) = A {0} (T fl (c)) 

3a' G R : T B (c) = T [a , } {x a {c)) = x a (T {a , +a} {c)) 

=>c = r {a , +a} (c), 
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and analogously 

3aeR:c = F {a} (c) => 3d G /? : c = A {fl / +fl} (c). 

Hence, by lemma [6j the cellular automata C and G are inverse to each other. ■ 
Thus only cell 0, the neighbors of cell 0, the cells in s$ and the neighbors of each cell in s# have to be 
considered when testing whether two fully asynchronous one-dimensional cellular automata are inverse 
to each other. 

Theorem 16. Phase space invertibility is decidable for fully asynchronous one-dimensional cellular 
automata. 

Proof. Let C = (X,N,Q,8) be an arbitrary fully asynchronous one-dimensional cellular automaton and 
T = {0} L)NL)£/ U (si + N). Consider each of the \Q\ 1^1 N fully asynchronous one-dimensional cellular 
automata G = (X,N,Q,y) in turn. Test whether for each pair c,c' G Q T of finite configurations with 
D c ,c' = {0} tne equivalence 

c ' = A {o}(c) <=^ c = r {0} (c') 
holds and for each finite configuration c G Q T the implications 

c = A{o} (c) =>■ 3a G 80 :c = Ts a \ (c) and 

c = r {0}(c) ^3aGi2/ :c = A {a} (c) 

hold. If this is the case for one automaton G, then C and G are inverse to each other, and C is invertible. 
Otherwise C is not invertible. 

Note that +N = {a + n \ a G s/,n G A^}. ■ 
Remark. It is in fact not necessary to consider all fully asynchronous one-dimensional cellular automata 
with the same neighborhood and states. One can show that if C is invertible, the local transition function 
of its inverse can be constructed easily from the one of C. Thus only one candidate needs to be considered. 

For further details see HU section 3.5]. 

6 Elementary Cellular Automata 

Using the decision algorithms it can easily be shown that exactly the purely asynchronous elementary 
cellular automata with Wolfram numbers 0, 35, 43, 49, 51, 59, 113, 115, 204, and 255 are invertible and 
that exactly the fully asynchronous elementary cellular automata with Wolfram numbers 33, 35, 38, 41, 
43,46,49,51,52, 54, 57,59, 60, 62,97, 99, 102, 105, 107, 108, 113, 115, 116, 118, 121, 123, 131, 139, 
145, 147, 150, 153, 155, 156, 195, 198, 201, 204, 209, and 211 are invertible. 
For further details see [ 1 1 , chapter 7]. 

7 Summary and Outlook 

We have introduced a definition of invertibility for asynchronous cellular automata, namely phase space 
invertibility, and shown that invertible purely asynchronous cellular automata are computationally uni- 
versal, that invertibility can be decided for arbitrary-dimensional purely and one-dimensional fully asyn- 
chronous cellular automata. 

It remains open whether invertible fully asynchronous cellular automata are computationally univer- 
sal and whether invertibility for higher-dimensional fully asynchronous cellular automata is decidable. 
If the latter can be shown to be true then invertibility of purely asynchronous cellular automata restricted 
to non-empty sets of active cells would also be decidable in higher dimensions. 
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